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PHYSICAL REVIEW B VOLUME 32, NUMBER 6 15 SEPTEMBER 1985 

Interaction between surface plasmons and localized plasmons 

G. S. Agarwal and S. Dutta Gupta 
* School of Physics, University of Hyderabad, Hyderabad-500 134, India 
(Received 8 March 1985) 

A theoretical analysis of the experiments of Holland and Hall [Phys. Rev. B 27, 7765 (1983)] is 
given. The metal-island film is characterized by an effective anisotropic dielectric function. Formu- 
las for reflection and transmission coefficients for an anisotropic layered structure are given. Calcu- 
lated normal and non-normal reflectivities are in reasonable agreement with experiments. The re- 
sults explain the observed shift in the surface-plasmon dispersion due to its interaction with local- 
ized plasmons. 



I. INTRODUCTION 

Recently Holland and Hall studied the interaction be- 
tween the surface plasmons in a metal film and the local- 
ized plasmons in a Ag-island film. 1 The shift in the 
surface-plasmon dispersion was measured as a function of 
the separation between the two metal films. 2 The charac- 
teristics of the localized plasmons 3 " 5 were studied by ex- 
amining the normal-incidence reflectivity as a function of 
the wavelength. The purpose of this paper is to show that 
the observed features can be explained in terms of a rather 
simple model in which the island film is described by an 
effective dielectric 5 function; though the results are very 
sensitive to the choice of parameters. Note that because 
of the structure of the islands, this film is anisotropic and 
thus in Sec. II we present the formulas for the calculation 
of the reflection and transmission coefficients from a lay- 
ered structure where each layer can be anisotropic. In 
Sec. HI we calculate numerically the reflection coefficient 
and obtain the changes in the dispersion characteristics of 
surface plasmons due to the interaction with localized 
plasmons. These results explain the broad features of the 
experiments of Holland and Hall. We also present ap- 



proximate expressions for the normal-incident reflectivi- 
ties which show explicitly the connection with the imagi- 
nary part of the effective dielectric function of the island 
film. 



H. REFLECTION FROM A MULTILAYER 
STRUCTURE WITH ANISQTROPY 

The experiments of Holland and Hall use a multilayer 
geometry consisting of a prism, a Ag film, a spacer layer 
of LiF, a. Ag-island film and air. The island film can be 
characterized by an effective dielectric function, which in 
general, is anisotropic. Hence we need to calculate the re- 
flection and transmission coefficients of a multilayer 
structure 6 with anisotropy. Such reflection coefficients 
can be calculated from the knowledge of the general struc- 
ture of the electromagnetic field in an anisotropic medi- 
um 
In 



In a uniaxial medium with dielectric tensor 
€ t =e xx =€ yy9 € n =€ a , cr a 0=O, a^£, and occupying the 
domain say L j <z <L 2f the fields are known 7 to have the 
general form 



c+vy) 



ut 



vt 



e, (u 2 +u 2 )' e t (w 2 +u 2 ) 



, 1 Ue to +U— -t,et-+*7) 



du dv , 



(2.1) 



(2.2) 



where 

S 2 = *i u Z_ v 2 t r 2 = ^. e> _iL (u 2 + „ 2) . 
C 2 C £ €n 

For each spatial component (i^u); + the general solution 
(2.1) involves four unknowns: e^, ej. The terms iej) 
correspond to ordinary (extraordinary) waves. Here s and 
t represent the z component of the propagation vectors for 
these waves. The boundary conditions at each interface 
involve the continuity of the tangential components of E 
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and H. It turns out to be convenient to match the com- 
ponents 

(uE x + vE y ), (vE x —uE y ), 

{vH x -UH y \ (uH x +vH y ) . 
We can now study the reflection from a layered struc- 
ture, using the usual matrix methods. 6 The matching of 
the field components E x , E y , H x and H y will yield four 
equations for four unknowns which in turn will result in 
4X4 matrices. It turns out that 4X4 matrices can be re- 
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duced to block-diagonal structure consisting of 2x2 ma- 
trices if we work with combinations 

{uE x + vE y ), [vE x ~uE y ) , 

(vH x -uH y ), -iuH x +vH y ) . 

From the structure (2.1) of the field, it is seen, for exam- 
ple, uE x + vEy—0 for the ordinary wave; whereas 
vE x — uE y ~0 for the extraordinary wave. Thus the ordi- 



nary and extraordinary waves get decoupled. 

Assume that zth uniaxial medium occupies the domain 
from z=L { _\ to L t . Matching the boundary conditions 
at z yields two equations 



fz = o or e 



(2.3) 



where m !. + ** 1 are 2 X 2 matrices: 



and 
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2_l„2 . r / e m 
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N 
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^0 



(2.6) 



Assuming incidence from the medium 0 (occupying 
domain 0 <z < — oo ), the reflection coefficient becomes 



(2.7) 



W d )l2 = 



(2.5) 



The reflection coefficient can be evaluated in terms of 
these m matrices. For a N layer system, we can write 
(2.3) as 



On the other hand, for incidence from the medium N + \ 
(occupying domain L N <z < oo ), the reflection coefficient 
is 



Assume that the Ag film occupies the domain 0<z <d u 
and that the spacer layer lies in d x <z <d 2 +d x . These 
two layers are characterized by the isotropic dielectric 
functions e x and e 2 . The island film occupies the region 
d x +d 2 <z <d x +d 2 -hd 3 . The prism is on top of the 
silver film. The dielectric function of the Ag film has 
been taken from the work of Johnson and Christy. 8 The 
dielectric function of the prism and LiF are chosen as 
6o=2.23 and e 2 = 1.9363, respectively. The dimensions of 
the metal islands are much smaller than the wavelength 
and their shapes are" generally ellipsoidal. Hence we use 
the effective-medium theory to write the approximate 
dielectric functions of the island film 



= 1 + 



q(e x -l) 



--t t n 



(3.1) 



(2.8) 



Thus given the direction of incidence and the dielectric 
properties of each layer, the reflection coefficients can be 
calculated. 

Note that the poles of (2.7) or (2.8) will determine the 
structure of the possible surface modes, in the layered 
structure, where each layer can be anisotropic. 

III. NUMERICAL RESULTS FOR REFLECTIVITY 
AND SHIFTS IN SURFACE PLASMON DISPERSION 

We now use the formulas of Sec. II to calculate the re- 
flectivity in the specific case used by Holland and Hall. 



where e x is to be taken from Ref. 8. 

The parameter q is the filling factor. The quantity / M 
depends on the geometry 4,5 of islands. It is well known 
for various shapes. However, the value of / to be used in 
(3.1) gets modified due to a number of other effects 3 such 
as due to the image dipoles, neighboring dipoles. Hence 
in what follows we leave f ^ as an open parameter and we 
try to estimate / M from the experimental curves of Hol- 
land and Hall. The effective thickness of the island film 
is given by 

d 3 = U m )/q , (3.2) 
where t m is the mass thickness. 
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FIG. I. The normal incidence reflectivity as a function of the 
wavelength of light incident from the island side of the struc- 
ture. The parameters chosen are e 0 = 2. 23, e 2 = 1 - 936, f, = 0. 1 9, 
q =0.04, -rfj = 50 nm, and the mass thickness of island film is 3 
nm. Different curves are for different values of the spacer layer 
thickness d 2 = (a) 40 nm (b) 55 nm, (c) 7*0 nm. For rf 2 = 55 nm, 
curve die) gives the reflectivity for mass thickness 1 nm (con- 
tinuous Ag film removed). 

We first examine the behavior of the normal incident 
reflectivity | R n | as a function of wavelength of the light 
incident from the metal-island side of the layered struc- 
ture. The reflectivity is calculated using Eq. (2.8) and by 
taking u = u — 0. From (2.2) it is clear that e t enters the 
calculations and hence the reflectivity will depend on e 3f . 
Thus the parameter /, will determine the normal incident 
reflectivity data. We have investigated | R„ | for a wide 
range of parameters. The results are shown in Fig. 1 . We 
find that the dip or the minimum on the blue side of 4000 
A is most pronounced for spacer-layer thickness ~ 55 nm. 
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FIG. 2. Real and imaginary parts of the effective dielectric 
function of the island film as a function of wavelength. The 
real and imaginary parts of the dielectric function of the con- 
tinuous Ag film (taken from Ref. 8) are also shown for compar- 
ison. 



Reasonable values of the parameters f t and q are obtained 
by fitting the normal-incidence reflectivity curves so that 
for a given q, the position of the resonance occurs at the 
reported value. The parameter q is chosen such that the 
optimum behavior with respect to d 2 is obtained. For 
these values of the parameters the effective dielectric 
function of the island film is shown in Fig. 2. The two 
maxima in Fig. 2 correspond to the two absorption mini- 
ma of Fig. 1. Figure 2 also shows for comparison the 
values of 6j obtained from the work of Johnson and 
Christy. Note that here we are considering a wavelength 
region in which Ag has a dielectric function whose real 
part is large and negative. In Fig. 1, we find that the ab- 
sorption minimum corresponding to. 4000 A is split into 
two minima. This happens because small values of q 
make d l large if the mass thickness is kept constant. We 
have verified that for smaller mass thickness such a split- 
ting is not there. This is also shown in Fig. 1 by the curve 
corresponding to mass thickness 1 nm. The extra feature 
is a consequence of the uniform-layer approximation for 
the island film. We now give an argument to explain the 
observed behavior in Fig. 1. We have to connect the ex- 
pression for the' reflectivity to the effective dielectric func- 
tion of the island film. The reflection amplitude, as ob- 
tained from (2.8), has the explicit form 



R = 



I- A 



A = 



m 22 + (€ 0 ) l/2 m 2} 



(3.3) 



with 



sin(fcicfi)cos(fc2<*2)sin(&3<*3) — 



— I sinUjd, )sin(/e 2 d 2 )cos(/c3cf 3 ) 

*2 



1/2 



cos(^i^i)sin(/c2^2)sin(^ 3 rf 3 )-fcos(/c 1 rf 1 )cos(k 2 d 2 )cos{k 3 d 2 ) , 



(3.4) 
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7 21 = - i —cos{k l d l )cos{k 2 d 2 )sin(kid 3 )--— ^cos(/c^, )sin(/c 2 <f 2 )cos(/c 3 ^3) 



+1 



1/2 



sinffcjrf, )sin(fc 2 */ 2 )sin(£ 3 ii 3 ) — 



(e,) 1/2 



sin(/: l c? 1 )cos(fc 2 ^ 2 )cos(/c3<i 3 ) 



m n = — i(€ ] ) l/2 sin{kid l )cos(k 2 d 2 )cos{k 3 d 3 )-{-i 



sin( /c 1 1 )sin( k 2 d 2 )sin( k 3 d 2 ) 



-/(e^^costfci^sinf^^ , 

1/2 

m n = cos( k 1 d 1 )cbs(k 2 d 2 )cos( fc 3 rf 3 ) — 



cos( ifc 1 1 )sin( & 2 <f 2 )sin( k 3 d 3 ) 
1 1/2 



sin( k 1 d x )sin( k 2 d 2 )cos( k 3 d 3 ) — 



The full expression is too complex. An understanding of 
the results of Fig. 1 can be obtained by analyzing the ap- 
proximate form. We make a series of approximations — 

(a) Ignore the presence of the spacer layer, i.e., set d 2 =0; 

(b) assume d 3 very small, so that we obtain reflectivity 
only to first order in d 3 ; (c) in the region of wavelength 
shown in Fig. 1 Ag has a large-negative real dielectric 
function, so we set {€\) xn ~in x \ (d) assume that the 
silver-layer thickness is large so that we can ^approximate 
sinU^cos^d,) by -<l/2i)e 1 \e l \ 
Under these approximations we can write 

[(€ o ) 1/2 -\]-ik 0 d 3 [e t -ie o ) 1/2 } " 
[(e o ) l/2 +l]-ik 0 d 3 le t + (e o )>"] 
for the case when the uniform silver film is absent and 



[eJ 1/2 -l]-^ 3 [6 r -( go ) 1/2 ] 
* 0== [(€„) 1/2 -f l]-//c 0 tf 3 [e, + (O 1/2 ] 



(3.10) 



if the silver film is present. On further simplifications, we 
find the following results for the reflectivity in the two 



1*0 



(6 0 ) ,/2 -l 



(e 0 ) ,/2 +l 



4k Q d 3 

;1 — lm£ r . 



(3.12) 



The approximate results [(3.11) and (3.12)] show the re- 
markable connection between the reflectivity and the 
imaginary part of the dielectric function of the island 
film. Thus all the resonances (localized plasmons) of the 
e t of the island film will be reflected as peaks in | R 0 \ 
and as dips in | R x | . This explains the numerical results 



sin(fcirfi )sin(fc 2 c? 2 )sin(A: 3 rf 3 ) , 



(3.5) 

(3.6) 

(3.7) 
(3.8) 



of Fig. 1. It should be borne in mind that the numerical 
results are based on the exact result (3.3). The change of 
sign in the two cases can be understood as in the wave- 
length region of interest, Ag has negative real dielectric 
function. Thus the Ag film acts like a mirror, which ac- 
counts for phase change on reflection. We have thus 
shown that the localized plasihon resonances appear as 
peaks in the light reflected from the island film. Howev- 
er, such resonances appear as absorption resonances in 
presence of the additional Ag film. 

Having obtained an estimate for q and we examine 
the coupling between surface plasmons and localized 
plasmons in the island film. Such a coupling is reflected 
in the reflectivity \R<f>\ data as a function of the angle of 
incidence </>. The light is considered to be incident from 
the prism side. The calculation is done using (2.7), with 
u =0 v =o/c sirup, and with incident light polarization 
such that ej inc = 0, e7 inc = 0, i.e., incident light is p polar- 
ized. Figure 3 shows | R# | as a function of <}>, both in 




tp (deg) 



FIG. 3. The reflectivity R + as a function of angle of in- 
cidence for light incident from the prism side. The parameters 
are /„=0.52, A.— 5145 A, and other parameters as in Fig. 1. (a) 
Ag, (b) Ag-LiF (55 nm), (c) Ag-LiF (55 nm) island. 
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FIG. 4. The dispersion curves k+ vs <j>. The solid (dashed) 
curves are in presence (absence) of the island film. Different 
curves are labeled by the thickness in nm of spacer layer. All 
other parameters are same as in Fig. 3. 



the absence and in the presence of the island film. The 
parameter /„ is estimated from the sum rule. The reflec- 
tivity minimum ^ is investigated for a large number of 
wavelengths and the results are plotted in Fig. 4. We thus 
find changes in the surface-plasmon dispersion relation 
due to coupling with localized plasma resonances. Figure 
4 gives X$ versus <f> for a range of the thickness of the 



spacer layer both in presence and;, absence of the island 
film. The shifts in the dispersion relation show an 'op- 
timum behavior with regard to the thickness of the spacer 
layer, which is in conformity with the experimental re- 
sults of Holland and Hall. It was also found that increas- 
ing q shifts the ^ versus <f> curves to the right. Similarly 
increasing anisotropy (i.e., decreasing /,) leads to the 
shifting of the dispersion curves to the left. 

IV. CONCLUSIONS 

In conclusion we have provided a theoretical model for 
the experimental observations of Holland and Hall. We 
have only tried to explain the broad features rather than 
detailed fitting with the data. The anisotropy of the is- 
land film is fully taken into account, though a simplified 
effective medium description of the island film is adopted: 
The model also assumes that the effects of substrate, 
neighboring islands can be accounted for by modifying 
the values of the parameters 9 f t and /„. A better 
knowledge of f n f„ and an improved description of the is- 
land film are expected to yield results much closer to ex- 
periments. 
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